Transport spectroscopy of Kondo quantum dots coupled by RKKY interaction 
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We develop the theory of conductance of a quantum dot which carries a spin and is coupled via 
RKKY interaction to another spin-carrying quantum dot. The found dependence of the differential 
conductance on bias and magnetic field at fixed RKKY interaction strength may allow one to 
distinguish between the possible ground states of the system. Transitions between the ground states 
are achieved by tuning the RKKY interaction, and the nature of these transitions can be extracted 
from the temperature dependence of the linear conductance. The feasibility of the corresponding 
measurements is evidenced by recent experiments by Craig et al^ 



The exchange interaction between a localized electron 
and itinerant electrons of a Fermi sea leads to the Kondo 
effect )2i Recently, the Kondo effect was observed in the 
quantum dot setting, where it causes an anomalously 
high conductance at low temperatures^ The itinerant 
electrons not only screen an impurity spin, leading to 
the Kondo effect, but also give rise to the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction between lo- 
calized spins.** The interplay between the Kondo screen- 
ing and RKKY interaction remains at focus of investi- 
gation of strongly correlated electron systems, and may 
play an important role in the heavy fermion metals.-^ This 
interplay is not trivial even in the minimal system allow- 
ing it, which consists of two localized spins "imbedded" 
into an electron Fermi seai^ 

Until recently, such two-spin system was the subject of 
theoretical investigations only. The quest for the prac- 
tical implementations of the quantum computing ideas 
has led to the interest in physics of spin devices, fn this 
context, transport properties of two spin-carrying quan- 
tum dots, coupled with each other by RKKY interac- 
tion, were studied experimentally* in a device schemat- 
ically shown in Fig. ^ The authors of Ref. Q were able 
to see the effect of RKKY interaction between the spins 
by monitoring the Kondo-enhanced conductance through 
one of the dots. Sufficiently strong RKKY interaction 
locks the localized spins into a singlet or triplet state and 
destroys or weakens the Kondo effect, thus liquidating 
the enhancement of the conductance. The experiments* 
demonstrated the fact of spin coupling and detailed quan- 
titative measurements seem to be within the reach of 
experimental capabilities. Further experiments will also 
bring an exciting prospect of an experimental investi- 
gation of the interplay between RKKY interaction and 
Kondo effect in a controllable setting. 

In this Letter we develop a theory of conductance spec- 
troscopy of the spin states of two s = 1/2 quantum dots 
coupled by the RKKY interaction. We start with a char- 
acterization of the ground state and low-energy excita- 
tion spectrum of the many-body system formed by the 
dots and itinerant electrons of a two-dimensional electron 
gas (2DEG). Next we concentrate on the case of a strong 
RKKY interaction, allowing us to treat the Kondo effect 




FIG. 1: The exchange interaction of S = 1/2 spins of dots 
Si and ^2 with the 2DEG results in the Kondo effect. The 
connection by weak contacts Ci and C2 to the bigger open 
dot M creates the exchange (RKKY) interaction between the 
two spins. The current J as a function of bias V is measured 
between source (s) and drain (d) contacts. 



perturbatively and elucidate the main features of the I- 
V characteristic of the device sketched in Fig. ^ The 
measurement of differential conductance G{V) = dl/dV 
enables one to distinguish between the singlet and triplet 
ground states of the two localized spins. We then con- 
sider the crossover between the singlet and triplet states 
on the linear conductance Go = G{V = 0). We identify 
signatures in the linear conductance which may allow to 
distinguish the singlet-triplet crossover from a quantum 
phase transition between these two states. Finally, we 
investigate the sensitivity of the RKKY interaction to an 
applied magnetic flux and determine the characteristic 
flux needed to change the sign of the RKKY coupling. 

The ground state of two localized spins interacting 
with each other and with a Fermi-sea of itinerant elec- 
trons depends sensitively on the relations between the 
corresponding interaction constants. Under very special 
conditions^ requiring fine-tuning of the system param- 
eters, a non-Fermi liquid state of the system may be 
reached. Away from these special points in the parameter 
space, the low-energy properties of the system are that 
of a Fermi liquid, and we will concentrate on this generic 
case. The ground state of the full system, including the 
localized and itinerant spins, is a singlet. Variation of 
the exchange couplings leads to a number of crossovers 
between different possible singlet states. The simplest 
Hamiltonian sufficient for describing the singlet states is 
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H = JiSiSi + J2S2S2 + J12S1S2 



(1) 



Here Si and ^2 are the spins of the two dots (iS'1.2 — 1/2), 
and Sa — J2(7<7' '4'a(j^<y<y''^aa' local spin densities 

of itinerant electrons. Electron operators V'qo- represent 
electron states^ coupled to a localized spin in dot a. 

First we note that at J12 = the Hamiltonian Eq. 
represents two independent s = 1/2 Kondo systems char- 
acterized by Kondo temperatures Tkq? oc exp(— 1/j/Jq.), 
where v is the band density of states at the Fermi level, 
a — 1,2. Atr = 0, each of the two local spins is 
fully screened by the corresponding modes of the itin- 
erant electrons. In these conditions, each of the dots 
provides a Kondo resonance for electron tunnelingi^ The 
independent screening, giving rise to resonances in tun- 
nelling through each of the dots, remains in effect for 
a sufficiently small inter-dot coupling. At stronger anti- 
ferromagnetic coupling, J12 ^ maxa {Txce}, the two local 
spins form a singlet of its own, and the Kondo resonances 
for the itinerant electrons vanish. In the case of ferromag- 
netic ( J12 < 0) coupling, the two localized spins form a 
triplet, which is fully screened by the itinerant electrons 
interacting with the two dots and the Kondo resonances 
for tunnelling through each of the dots persist. 

The above consideration shows that the zero- 
temperature linear conductance changes from a large 
value (induced by the Kondo resonance) at (— J12) ^ 
maxalTxa}, to a small value at J12 ^ maxclTKa}. In 
the over-simplified representation of the device sketched 
in Fig. n by the Hamiltonian J^l, these asymptotes are 
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and G = 0, respectively (here Gs and Gd are the 
conductances of the junctions connecting the dot with 
the source and drain leads, respectively). In the spe- 
cial case^ of Tk2 = 0, the transition between the two 
asymptotes is a jump at J12 — between two Fermi 
liquid states. If Tk2 is finite, the transition shifts to 
J12 ~ Jki/ln(rKi/Tk2) for Tk2 < Tki, and remains 
at positive J12 ^ maxalTKa}. This transition occurs via 
passing through a non-Fermi liquid state, which belongs 
to the same universality class as the two-channel s = 1/2 
Kondo problemiS. The existence of such non-Fermi liq- 
uid state hinges on a special particle-hole symmetry'' of 
the Hamiltonian However, in the case of generic pa- 
rameters of a quantum dot the Hamiltonian of the sys- 
tem also includes other terms {e.g. potential scattering 
terms leading to the elastic co-tunnelling) which violate 
the required symmetry. In this case the quantum phase 
transition between the two Fermi liquids is replaced by a 
smooth crossover The zero-temperature conduc- 
tance varies smoothly and monotonically with J12, as 
shown by the solid line in Fig. [5fa).^^ 

In the absence of the Zeeman splitting, the energy of 
exchange interaction J12 sets the threshold for the in- 
elastic electron scattering, accompanied by a change of 
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FIG. 2: (a) Linear conductance Go as a function of the RKKY 
coupling constant J12 at T = 0. The crossover from the 
triplet state J12 < at |Ji2| ^ maxa{rKa} to a singlet state 
J12 ^ maxa{rKa} takes place at |Ji2| maxa{rKa} and 
results in the decrease of Go from a value close to the unitary 
limit Gu to a much smaller value controlled by the elastic co- 
tunnelling. (b,c) Differential conductance G{V) at T = and 
strong RKKY coupling, |Ji2| 3> maxalTka}. Solid lines: the 
second-order in exchange Ji result; dashed lines: schematic 
representation of the enhancement of G{V) resulting from 
higher-order in Ji calculation. 



the spin state of two dots. Far above the threshold, 
|eT^| I •^12 1, the processes with flip of the spin are 
allowed and, in the leading order, the differential conduc- 
tance G{y ) coincides with the conductance of a single- 
quantum-dot device: 



G{V) 
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Here factor 3/4 corresponds to the square of the operator 
of spin S ~ 1/2, and the logarithmic term represents the 
Kondo renormalization of the exchange interaction in the 
weak coupling limit, eV ^ Tki- 

Conductance Eq. (jSjl arises from the scattering ampli- 
tude evaluated to the lowest order perturbation theory in 
the renormalized interaction constant. Within this accu- 
racy, conductance has a sharp step at \eV\ = \ Ji2\- The 
height of the step depends on the sign of J12, i.e., on 
wether the coupling of the two localized spins is ferro- 
magnetic or antiferromagnetic. In the case of ferromag- 
netic couphng, J12 < 0, conductance G[V) is reduced by 
factor 3/2 from the value given by Eq. Q if the bias is 
lowered just below the threshold \eV\ — —J12. The renor- 
malization group analysis of the Hamiltonian |^ shows 
that after the reduction at \eV\ = — J12, the differen- 
tial conductance again monotonically increases as bias V 
decreases: 



GiV) = 
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W \eV/TKt\ 
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Here factor 2 corresponds to the square of S" = 1 spin 
operator and the logarithmic increase of G{V) reflects 
the S = 1 Kondo effect with the triplet Kondo tempera- 
ture Txt. In case rK2 *C Iki <^ (— J12) an evaluation of 
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FIG. 3: The contour plot of fi{Ji2,V, B) is shown in the 
bias (V) - magnetic field (B) plane. Function /i determines 
the difFerential conductance G{V) at strong RKKY coupling 
|Ji2| > ma.Xa{TKa}, see Eq. @. 



the first logarithmic correction to the conductance gives 
Txt = T^i/\Ji2\, i-e. at large |Ji2| the triplet Kondo 
temperature is much lower than uiaxaiTKa} ■ For the an- 
tiferromagnetic coupling between the two localized spins 
the particle-hole symmetric Hamiltonian yields zero 
conductanceii^ 

Magnetic field may shift and split the step in conduc- 
tance G{V). The proper generalization of Eq. Q to in- 
clude the Zeeman splitting gfiB by magnetic field B reads 



G{V) 
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with a step-like function /i, presented in Fig. 01 At 
gfiB :s> I J12I, \ eV\, the locahzed spins are aligned along 
the field and /i = 1/4. If J12 < 0, Eq. (0) must be 
supplemented by condition ma,x{g ^B , eV} }^ \/T^^Jv2\ 
otherwise a S* = 1 Kondo effect develops, see Eq. Q). 

To assess the applicability of the above results ob- 
tained in the Born approximation for the electron am- 
plitude of tunnelling through the dot, one may evaluate 
the next-order correction to that amplitude. This cor- 
rection diverges logarithmically at \eV\ = \Ji2\, similar 
to the divergence occurring in the scattering off a sin- 
gle localized spin in the presence of Zeeman splittingii^ 
The divergence is cut off by cither the temperature!^ or 
the relaxation rate of the localized spins in the out-of- 
equilibrium conditions^ As the result, steps in the dif- 
ferential conductance, obtained within the Born approx- 
imation, become asymmetric maxima of a finite width. 
Here we present an estimate for the correction to Eq. Q 
only for the case of strongly asymmetric setup, e.g., 
Gs/Gd 1, and assuming the lowest temperature, which 
allows for the sharpest steps in the differential conduc- 
tance. In these conditions, the spin relaxation ratel- re- 
sponsible for the cut-off of the divergency is of the order 
of 1/t, - |Ji2|/[ln^(|Ji2|/TKi)]. The widths of the fea- 
tures replacing steps in the differential conductance are 
of the order of ~ 1 /t^ , and the amplitudes of the correc- 
tions to steps described by Eq. lO are 
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The amplitude of the correction is small if | J12I ^ Tki, 
and the main effect of the higher-order terms of the per- 
turbation theory is the smearing of the steps. In order 
to determine the nature of the ground state (singlet or 
triplet) by conductance measurements in the presence of 
magnetic field, the Zeeman splitting energy must exceed 
1/ts of the features in the differential conductance. The 
overall bias dependence of the differential conductance at 
_B = is sketched in Fig. I3b,c). 

Linear conductance Gq of the dot Si is determined by 
the scattering T-matrix Ti (e) for e ^T: 



Go(T, J12) = Gu J 



ImTi{e)de 
4Tcosh2(e/2r)' 



(7) 



In the unitary limit Ti(0) = i. For e ^ T ^ maxalTxa} 
we may use the Born approximation^^: 

Imri(e) = v^J^ J K{e - e')l±^de', (8) 

where K{uj) = 2ttJ2^^, p^\{C\Si\0\^ 6{iu + - E^,) is 
calculated with respect to the exact quantum states |^) 
of the system of two localized spins, E^ is the energy of 
state 1^), and is the density matrix cx exp{—E^/T). 
The exchange constant J7e is logarithmically renormal- 
ized by the Kondo effect: = l/[:y\nE/TKi] with 
S = max{£,T, IJ12I}. 

Within Born approximation, the linear conductance 
depends on J12 only through the ratio J12/T, 
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ln^(T/TKi) 
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The dependence of Go(T, J12) on J12 has a maximum 
in the region |Ji2| ^ T and conductance is higher on 
the triplet side of the crossover. (Note however, that 
the exponentially small value of conductance far in the 
singlet region is an artefact of our model.J'-) At nega- 
tive J12, as temperature T decreases, the conductance 
Go(T, J12) grows and reaches Gu at T < T^^/lJu]. The 
shape of the Go(T, J12) vs. J12 eventually approaches 
the step sketched in Fig. In the generic case of the 
zero-temperature crossover between the two fermi-liquid 
states, the width of the step saturates and remains finite 
in the limit T — )■ 0. This width is not universal and de- 
pends on the terms in the exact Hamiltonian beyond the 
approximation of Eq. If, however, the parameters 
of the system are tuned properly, and the variation of 
J12 takes the system through the non-Fermi-liquid state 
at certain value of J12, then the step width decreases 
with lowering the temperature as Vt, see, e.g., |2Clj |. Fi- 
nally, the limit J2 = within the model Eq. corre- 
sponds to a sharp transition between two Fermi-liquid 
states at T = 0. A straightforward analysis, similar to 
Ref. 1^, yields the estimate Tki / ln(TKi/r) for the tran- 
sition width at finite temperature. 

The contact interaction of the localized spins with the 
itinerant electrons of dot M at points Gi and G2, see 
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Fig. n results in the indirect exchange interaction be- 
tween the locahzed spins. Unhke the textbook RKKY 
interaction facilitated by freely propagating electrons^, 
here magnitude and sign of J12 are random, reflecting the 
chaotic electron motion in the dot M. Roughly, the typ- 
ical value of J12 is ^ SiJciJc2 7 and the typical magnetic 
flux <i>c needed for changing J12 substantially is ~ $0. 
Here Jci 2 the diniensionless constants of the contact 
exchange interaction at points Ci^2', <5i — (I'A)^^ is the 
mean level spacing of one-electron energy levels in dot 
M, I' is density of states of the 2DEG, A is the area of 
the dot M, and $0 = hc/e. 

The RKKY coupling J12 may be expressed in terms 
of the scattering matrix 5i;2 of an electron propagating 
from contact Ci to contact C2: 

f de 

J12 = - 2JcM — n(e)Im{52;i(e)5i;2(£)} . (10) 

Here n(e) is the Fermi function and we assumed that elec- 
tron propagation in M is spin independent: Sia-2(7'{£) = 
(5o-o-'5i;2(e). Within the random matrix theory, it is re- 
lated in a standard way^-'^ to the one-electron Hamilto- 
nian 7i of dot M . The ensemble average ( J12) — 0, and 
we calculate the correlation functionSS of the RKKY con- 
stant Ji2{^) over realizations of matrix 'H{^) at two val- 
ues of magnetic flux $1^2 threading dot M, 

(Jl2(*l)Jl2($2)) = ^2AJc^^ • (11) 

Here £± — 7csc + ± <i>2)^/<&§; numerical factor 



K '--^ 1 depends on geometryjSLS^ ^^^^ = N5i/(2n) is 
the electron escape rate from the middle dot into 2DEG 
through N open channels; i?Th = vf/"/A is the Thou- 
less energy; and wp is the Fermi velocity. At T > f± 
should be replaced by T . 

Equation shows that the RKKY constant is sym- 
metric with respect to the inversion of magnetic flux 
<f> — <i>, and yields the correlation flux value '^^^'^ ~ $0. 
The flux ^^^"^ is much larger than the correlation flux 
<f>^ ^ ^7osc/-£'Th'&o for the conductance of an open 
quantum do^. The difference between <^'l^^ and $^ oc- 
curs because the contribution to J12 originates from en- 
ergy levels within the spectrum "window" E'tii, whereas 
the conductance of a dot is usually determined by levels 
within much shorter energy interval 7esc- 

In conclusion, the presented results may help one to 
determine the spin states of quantum dots coupled by 
RKKY interaction from transport measurements, see 
Eq. ^ and Fig. The evolution of the linear con- 
ductance with the variation of RKKY interaction con- 
stant allows one to follow the transitions between vari- 
ous ground states of the system. Finally, we found that 
the magnetic field flux needed for variation of the RKKY 
coupling is of the order of flux quantum $o- 
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